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Abstract 

We have utilized the concept of valon model to calculate the spin structure func- 
tions of proton, neutron, and deuteron. The valon structure itself is universal and 
arises from the perturbative dressing of the valence quark in QCD. Our results agree 
rather well with all of the relevant experimental data on g^'^^'^ and and suggests 
that the sea quark contribution to the spin of proton is consistent with zero. It also 
reveals that while the total quark contribution to the spin of a valon, AS^aion, is 
almost constant at Q'^ > 1 the gluon contribution grows with the increase of and 
hence requiring a sizable negative orbital angular momentum component L^. This 
component along with the singlet and non-singlet parts are calculated in the Next- 
to-Leading order in QCD . We speculate that gluon contribution to the spin content 
of the proton is about 60% for all values. Finally, we show that the size of gluon 
polarization and hence, Lz, is sensitive to the initial scale Qq. 

1 INTRODUCTION 

A central goal in the study of QCD is to understand the structure of hadrons in terms of 
their quark and gluon degrees of freedom. The most direct tool and sensitive test for prob- 
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ing the quark and gluon substructure of hadrons is the polarized Deep Inelastic Scattering 
(DlS)processes. In such experiments detailed information can be extracted on the shape 
and magnitude of the spin dependent parton distributions, Sqf{x,Q'^). Deep inelastic 
scattering reveals that the nucleon is a rather complicated object consisting of an infinite 
number of quarks, anti-quarks, and gluon. It is a common belief that other strongly inter- 
acting particles also exhibit similar internal structure. However, under certain conditions, 
hadrons behave as if they were composed of three (or two) constituents. Examples are 
the magnetic moments of the baryons, meson and baryon spectroscopy, the meson-baryon 
couplings and the ratio of total cross sections such as ^(^'^^^ and so on. Thus, it seems to 
make sense to decompose a nucleon into three constituent quarks called U and D. They 
would carry the internal quantum numbers of the nucleon. On the other hand, in DIS 
one observes that a nucleon has a composition of essentially an infinite number of quark- 
antiquark pairs and gluons, in addition to its valence quarks. One might identify the 
valence quark with a constituent quark, but this would imply that the three quark picture 
is a very rough approximation and both qq pair and gluon degrees of freedom need to be 
added to the picture. In doing so, it would be very difficult to understand why the three 
quark picture of a baryon works so well in many circumstances. One way of reconciling 
this apparent contradiction is to consider a constituent quark as quasi-particle with a non- 
trivial internal structure of its own; consisting of a valence quark and a sea of qq pairs and 
gluons. Such an interpretation of constituent quark is not new, more that 30 years ago it 
was advocated by Altarelli and Cabibo [1]. R.C. Hwa developed a more elaborated version 
by introducing the so called valon model [2] (the term which we will use hereafter) and 
applied it to a variety of phenomena with great success. More recent indication for the 
existence of the valon can be inferred from the measurements of the Natchmann moments 
of the proton structure functions at Jefferson Laboratory. They point to the existence 
of a new scaling that can be interpreted as a constituent form factor consistent with the 
elastic nucleon data [3]. This finding suggests that the proton structure originates from 

2 



elastic coupling with extended objects inside the proton. In References [4] [5] [6], the valon 
concept is utilized to calculate the unpolarized structure function of a number of hadrons. 
The results are in excellent agreement with the experimental data. Altarelli has also cal- 
culated the pion structure function in the constituent quark representation [7], using a 
deconvolution procedure. On a more theoretical front, M. Lavelle and D. McMullan [8] [9] 
proved that one can dress a QCD Lagrangian field to all orders in perturbation theory and 
construct a constituent quark in conformity with the color confinement. Prom this point 
of view, a valon is defined as a structureful object emerging from the dressing of a valence 
quark with gluons and qq pairs in QCD. Chiral Models in the realm of non-perturbative 
QCD also require the dressing of a valence quark and thus producing a structureful object. 
These results and the success of the valon model in describing the unpolarized structure 
of hadrons and a number of other low Pt hadronic phenomena lend credit for the study 
of hadrons in the valon representation. 

In this paper we calculate the polarized structure of a valon and extract the polarized 
structure function, gi, of the nucleon and the deuteron and compare the results with the 
experimental data. A number of similar attempts are also made to derive the nucleon 
spin from the quark models [10], [11], [12], [13], [14]. Our model differs from those in that 
we calculate the polarized structure of a constituent quark (the valon) directly from QCD 
processes in the Next-to-leading and investigate its peculiarities and distinctive features. 
In order to be clear, we define a valon as a valence quark plus its associated sea partons, 
emerging from dressing processes. In a bound state problem those processes are virtual 
and a good approximation for the problem is to consider a valon as one integral unit whose 
internal structure cannot be resolved. Therefore, it is assumed that the spin of the nucleon 
is provided by the combination of the spins of the valons. In a scattering situation, on 
the other hand, the virtual partons inside a valon can be excited and be put on mass 
shell. It is therefore more appropriate to think of a valon as a cluster of partons with some 
momentum and helicity distribution. If the valon has a non-trivial internal structure then 
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the question arises whether its spin structure is also complex as it seems to be the case 
for the nucleon. This issue will also be addressed. In Reference [15] this model is used 
to calculate polarized structure functions. While obtaining results that are in agreement 
with the experimental data, however, it contains misleading and at points even counter 

intuitive ingredients. We will address them throughout this paper. 

The organization of the paper is as follows: First we will outline the formalism for calcu- 
lating the spin structure of the valon, then the polarized structure function of the nucleon 
and deuteron will be evaluated and it will be shown that the orbital angular momentum 
of partons in a valon plays a central role in describing the spin of nucleon. 

2 Polarized Valon Structure 

In this section we will utilize the extended work done on the development of NLO calcula- 
tion of the moments, to evaluate the polarized structure of a valon. We should stress that 
this is not a new next-to-leading order calculation, but it is an exploration of the existing 
calculations in the valon framework. 

By definition, a valon is a universal building block for every hadron; that is, its structure 
is independent of the hosting hadron. The valons play a dual role in hadrons: (i) they 
interact with each other in a way that is characterized by the valon wave function and 
(ii) they respond independently in an inclusive hard collision with a dependence that 
can be calculated in QCD at high Q^. In role (i) they are the constituents of bound state 
problem involving the confinement at large distances. In role (ii) they are quasi-particles 
whose internal structure are probed with high resolution and are related to the short dis- 
tance problem of current operators. This picture suggests that the structure function of a 
hadron involves a convolution of two distributions: valon distribution in the hadron and 
the parton distribution in the valon. In an unpolarized situation we may write: 




(1) 



4 



1 J X y y 



valon 



where F2"^'''"(|, Q^) is the structure function of the probed valon and can be calculated 
in Perturbative QCD to a certain degree of approximation. If is small enough we may 
identify F'"°-''°'^{x,Q'^) as 5{z — 1) at some point, for the reason that we cannot resolve its 
internal structure at that value. Similarly, for a polarized hadron we can write 

where SG'^^i^^{y) is the helicity distribution of the valon in the hosting hadron and 
gvalon ^ q2'^ is the polarized structure function of the valon. At high for a U-type 
valon one can write g'"°'^°'^ as follows: 

"^91 (z, Q^) = U^Gu + 5Gu) + ]-{6G d+5Ga+ SGs. +6G±) + ... (3) 

where all the functions on the right-hand side are the helicity functions for finding polarized 
quarks with momentum fraction z in a U-type valon at that Q^. These functions or certain 
combinations of them can be calculated in QCD and assumed to be known as we will meet 
them later. Similar expression can also be written for the D-type valon. To describe the 
polarized parton distribution inside a valon , we will work in the moment space, where 
the moment of the polarized parton distribution in a valon is defined as: 

A/(n,Q2)= j\^-^5f{z,Q^)dz. (4) 
Jo 

Sf{z, Q^) corresponds to parton helicity densities in a valon. The moments of the valon 
structure function are expressed completely in terms of through the evolution param- 
eter t: 

t = ln^. (5) 

We work in MS scheme where ^qcd given by 

A& = .«xp{-i[^-|M^ + |)l}. (6) 
where ji is the factorization scale and the (3 functions are as follows 

/3o = y - /^Ti., p^ = ^cl-f'^CATF-^fCFTp. (7) 



Here Ca = 3, Tj? = ^, Cp = |, and / is the number of active flavors. A NLO fit to the 
gi/Fi with massless quarks is performed in [16] and favors a value Aqcd = 0.235 ± 0.035 
GeV. This is very close to our choice of A.qcd = 0.22 GeV for the unpolarized case, we 
will maintain this value along with Qq = 0.283 GeV'^ as in [4]. 
Moments of the polarized valence and sea quarks in a polarized valon are: 

5M s,. =6MNs(n,Q^) (8) 

valon 

SMs^ = ^{SMs - SMNs){n, Q^) (9) 

where / is the number of flavors and 5Ms,ns are polarized singlet and non-singlet moments 
deflned as: 

5M,s.in,Q^) = {1 + "-(Q^) -"-(Q^ (_g)(,pa)n^ _ A<)")}L-^i)^-^'«°" (10) 



^ 27r 2tt ^ 
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The column matrix on the right hand side describes our initial input densities and they 
constitute an essential part of this work. We have inferred them from the properties of 
the model: The valon structure function has the property that it becomes S{z — 1) as 
is extrapolated to Qq (beyond the region of validity). This mathematical boundary 
condition means that the internal structure of the valon cannot be resolved at Qq in 
the NLO approximation. Consequently, when this property is applied to Eq. (2), the 
structure function of the nucleon becomes directly related to xSC'^^i^^^ (x) at that value of 
Qq] that is, Qq is the leading-order effective value at which the hadron can be regarded 
as consisting of only three (two) valons for baryons(mesons). In the moment space it is 
the Mellin transform of the (5-function, being equal to one, that enters. Naturally, then 
SMg{n, Qo) = and it is reasonable to set 5f{z, Qq) = f{z, Qq), for the quark sector. As 
Q^ increases beyond a small enough value, say Q^, where we may identify valon structure 



as S{z — 1), one expects the valon structure to develop a tail in the < 2; < 1 region 
due to gluon radiation. Reliable calculations are only possible for higher Q"^ values; let 
it be for Q'^ > Q\. Between Ql and Q\ higher twist terms are involved and the picture 
is more complicated, but that is also the region where the most important part of 
evolution takes place. In reference [15] the initial moments of both gluon and singlet sectors 
are set equal to one. Such a choice has no justification in the valon model. Moreover, 
it is known that a fully saturated initial input density for the gluon is disfavored [16]. 
Basically, initial input densities are determined from the experimental data, but for the 
valon there is no experimental data, therefore our choice of initial input densities, based 
on the mathematical conditions of the model does not invalidate it, nor does it violate the 
positivity constraint. 

In Eqs.(10,ll), L = as((5^)/as(Qo)' ^'^'^ JpW'^ is 2x2 singlet matrix of splitting functions, 
given by 



^p(0)n 



(12) 



where SPj-^^ are the n*'* moments of the polarized splitting functions and U accounts 
for the 2-loop contributions as an extension to the leading order. The explicit forms of 
these functions are given in [17] in the next-to-leading order. Now it is straightforward 
to calculate the moments of polarized partons inside a valon at any or t value. These 
moments are shown in Figure 1. The z-dependence of the polarized parton distributions 
is obtained by utilizing the usual inverse Mellin transformation. 

SQnss ci^, Q^) = - Im[e''^z-^-'"^'^6M^^'^'^{n = c + we''^, Q^)]dw, (13) 

TT Jo 

where NS, S, and G stand for non-singlet, singlet and gluon, respectively. In what follows 
we shall only be interested in quantities averaged over z and thus in the n = 1 moments. 
The first moments are defined by 

A/(Q')= / dzSf{z,Q^) (14) 



Figure 1: Moments of singlet, non-siglet quark, and gluon distributions in a valon at several values. 

There is a simple physical interpretation for the quantities like AS, and : they 
are related to the total z component of quark and gluon spins, thus 

< Sz >q= ^M, < S, >g= Ag (15) 

These quantities for the valon are shown in Figure 2. The results imply that the total 
quark contribution to the spin of a valon, decreases from 1 to AS = 0.88, in the range of 
= [0.283, 1] GeV^ and remains almost independent of thereafter. In other words, 
if the valon consisted only of quarks (valence plus sea) it would have been enough to 
account for 2± 90% of the valon spin at > 1. Evidently, however, there is a sizable 
gluon component, Ag, which increases with as shown in Figure 2. In figure 3 we have 
shown the variation of Aqgea and Aqyaience 

as a function of Q^. The variation of Aqg^a 
and AS with is very marginal, we have checked that Aqyaience = [1; 1-08] for the range 
of Q2 = [0.283, 10^] GeF^; whereas Ag 

sea varies from to -0.043 for the same range of 
Q^. This weakly dependent behavior of Aq^aience and Aqsea, are well understood: in 
the leading order one expects the total quark spin to be constant due to the vanishing of 
quark anomalous dimensions at n = 1. In the Next-to-Leading Order, however, they are 
marginally dependent due to APj^g ^ 0. The fact that sea quark polarization in the 
valon is consistent with zero can also be understood on theoretical grounds. The valon 
structure is generated by perturbative dressing in QCD. In such processes with massless 
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Figure 3: variations of t^qvaimce. and Agsea in a valon as a function of . Pull squares are the results 
calculated with Q\ = 0.283 GeV^ and open circles correspond to Qo = 1 GeV^ . 

quarks, helicity is conserved and therefore, the hard gluons cannot induce sea polarization 
perturbatively. It is also worth to note that /^qsea — is in good agreement with HERMES 
data [18] [19]. 

Had we chosen Q'q = 1 GeV"^, we would have obtained AS = 1 for all values of Q"^ 
with little change in Aqsea, and with a much reduced gluon polarization, as compared to 
our results with Q§ = 0.238 GeV"^ as can be seen in Figure 3. Such a choice with initial 
gluon helicity distribution equal to zero, however, is inconsistent with the mathematical 
condition of the model. The use of a different initial gluon helicity distribution, instead 
of zero, also would have been a pure guess and at best would require data fittings that we 
have avoided. 

It is obvious that due to large gluon polarization in a valon, these results do not add up 
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to give the spin ^ of the valon and do not satisfy the sum rule 5 = 5 AS + Ag. The gluon 
contribution to the valon spin grows as increases, while AE remains almost unchanged 
beyond = 1 GeV^. It is the present wisdom that the above sum rule must be replaced 
with a more realistic one 

1 = ^AS + A5 + L, (16) 

where, is the orbital angular momentum carried by the sea partons {q — q pairs and 
gluons) within the valon. The size of this orbital angular momentum turns out to be 
large and negative, mainly competing with the gluon contribution. Ratcliffe [20] was the 
first to point out the necessity of including orbital angular momentum dependence of the 
evolution equation and predicted a negative value for < L >2 of the sea partons in the 
proton. In Fig. 4 we present the Next-to-Leading Order result for the orbital angular mo- 
mentum in a valon, L^"'°"((5^). A leading order calculation is given in [21]. The existence 
of non-zero internal orbital momentum in the valon implies that there are substantial 
correlations among partons. It has been argued that the presence of quark pairs inside 
hadron resembles superconductivity [22] . An extension of the theory of superconductivity 
to the anisotropic case [23], shows that the presence of anisotropy leads to axial symme- 
try of pairing correlations around the anisotropy direction and to the particle currents 
induced by the pairing correlations. The particle number conservation, then requires that 
the cloud of correlated particles to rotate around the central particle in opposite direc- 
tion, the so-called Backflow current. This pairing induced orbital angular momentum is 
proportional to the density of correlated particles. The analogy seems to match the valon 
picture, where the internal structure of the valon originates entirely from QCD processes. 

3 Polarized Nucleon Structure Function 

In the previous section the polarized structure of a valon is completely specified. We are 
now in a position to carry forward and investigate the implications of the model at the 
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Figure 4: Orbital angular niomcnturn, L""'°"(Q^), component of partons in a valori as a function of . 
The results corresponds to Ql = 0.283 GeV'^. 
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hadronic level. Our starting point is Eq.(2), where the only unknown element is the valon 
helicity distribution, (^G^^;^^ (y) , since the valon structure function ^^"^o" is now completely 
given. In the analysis of [21], for the leading order, we assumed that the polarized valon 
distribution is related to the, by now well determined, unpolarized valon distribution via: 

5G^{y) = 5F^{y)G^{y) (17) 

where, Gj{y) is the unpolarized valon distribution of j = U,D kinds. Gj{y) are given in 
[4] [5] [6] for a variety of hadrons. They mimic the hadron wave functions. 
In the absence of experimental knowledge on <^G^^;Q„(y), the safe way to determine them 
is to fit the experimental values of at some with the form given by 

5Fj {y, Ql) = Njy'^i (1 - yfi (1 + a.yO-^ + hjy + cjy^-^ + djy'') (18) 

where, j stands for U and D type valon and parameters aj, (5j, etc. are given in Table I. 
Consequently, the polarized valon distributions in a proton are now completely specified 
and are given by: 

8Gu{y) = dFuGu SGoiy) = SFdGd (19) 



Table I. Numerical values of the parameters in Eq. (18). 



valon{j) 


Nj 


aj 


Pj 


aj 


bj 




dj 


U 


3.44 


0.33 


3.58 


-2.47 


5.07 


-1.859 


2.780 


D 


-0.568 


-0.374 


4.142 


-2.844 


11.695 


-10.096 


14.47 



Now we can calculate the polarized hadronic structure function, gi- For this purpose 
all that we need is to substitute SGu{y) and SGoiy) from Eq. (19) into Eq. (2) and 
perform the convolution integral. In Fig. 5 we present the results for proton, and 

compare them with the experimental data [18], [19] [24], [25], and with the calculations of 
[16], [26], and [27]. 
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Figure 6: x5uv{x) and xSdv{x) at — 2.5 GeV^. The curves are the model results and the data points 
are from Ref. [18, 19] 

It is evident that the model calculation is in good agreement with the experimental 
data. In Fig. 6 we compare the calculated results of x5uy{x) and x6dy{x) with the 
experimental values from HERMES at = 2.5 GeV^. Results from other analysis are 
also shown. 

As a further comparison, in Figs. 7 and 8 we present results for xg'^ and xgf. These 
results are also in agreement with the experimental data [24] and the analysis of References 
[16], [26], and [27]. 



4 First Moments and the Spin of Proton 

The first moment of polarized proton structure function, defined by 

r?= C cf,{x,Q^)dx (20) 

^0 
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Figure 7: Polarized neutron structure function, xg^, as a function of x at = 2,3, and 5 GeV^. Data 
points are from [24]. 




Figure 8: Polarized deuteron structure function, xgf, as a function of a; at = 2, 3, and 10 GeV^. Data 
points are from [24]. 
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can be related to the combinations of the quark spin components via 

r? = ^ E = IE4<P^'\ «7M75g I P, ^ > s^- (21) 

q q 

Our results for are listed in Table II. The moments of the polarized quark contributions 
to the spin of proton at, say, = 3 GeV^ are 

Au^alence = 0.820, Ad^alence = -0A22, Ag,,„ ~ (22) 

For other values of similar results are also obtained: for = [2,10], we have 
Auy = [0.816,0.827] and Ad^ = [-0.420,-0.426]]. These results for the first moment 
of the quantity r^'"'*^, yield the values that are presented in Table II. 



Table II. Numerical values for at several values 



-pN 
1 


g2 = 2GeF2 


g2 = 2.5GeF2 


Q2 = 2,GeV^ 


Q2 = ^QQy2 


Q2 = lOGeV^ 


P 


0.1132 


0.1153 


0.1168 


0.1200 


0.1245 


n 


-0.0630 


-0.0459 


-0.0554 


-0.0548 


-0.0546 


d 


0.0305 


0.0386 


0.0341 


0.0357 


0.0377 



The available experimental values for F^'" are obtained in a range of Q^, rather than 
at a fixed and a fit to all data at = 5 yields Ff = 0.118 ± 0.004 ± 0.007 and 
= -0.048 ± 0.005 ± 0.005 [28]. More recent data from HERMES [31] suggest that 
Ff = 0.1211 ± 0.005 ± 0.008 and Ff = 0.0436 ± 0.0012 ± 0.0018 at = 5 GeV2. The 

values given in Table II matches these experimental results. 

Additional information on the quark polarization is also available from the low-energy 
nucleon axial coupling constants g\ and g\: 

g\ =<p,s\ u-fi^'y5U - d-fi^j^d \p,s> s'' = Au{Q'^) - Ad^Q"^), 

_ (23 
g\ =<p,s\ U7^75« + d7^75ci - 2s7/.75S \p,s>s^' = Au{Q^) + Ad{Q^) - 2As{Q'^). 

Since there is no anomalous dimension associated with the axial-vector currents, and 
A^, the non-singlet couplings, g\ and g^ do not evolve with and hence can be de- 
termined from low-energy neutron and hyperon /3-decays. The experimental values are 
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g\ = 1.2573 ± 0.0028 and g\ = 0.579 ± 0.025. The constraining values of gY are the 
ones used in most analysis performed in order to fix sea quark contributions, in particular 
that of AS'. We have not considered this restriction a priori; instead we want to see if 
the results of the model can reproduce these numbers. Prom the stated values for Aw^ 
and Ad^ we obtain g\ = 1.240 — 1.253 which accommodates the experimental value with 
an accuracy of 2%. Our agreement with g\ is not as good as g\. However, this does not 
invalidate the results of the model, for there are serious objections [29] [30] ( mainly due 
to mu,d <.rns) to 

g\ = Au + Au + Ad + Ad- 2{As + As) = 3F - D = 0.579 ± 0.025, (24) 

in contrast to the unquestioned isospin SU{2) symmetry (m^ ~ m^) that gives rise to 
the value of g\. Experimental data from HERMES also puts the value of g^ at 0.274 it 
0.026 ± 0.011 in the range of 0.02 < a; < 0.6 [31] which is substantially less than the value 
inferred from hyperon decay. We have obtained g\ = 0.39. Findings of [31] also suggests 
that As + Ss is consistent with zero. The HERMES data [31] in the measured region of 
X > 0.02 gives As + 6s = 0.006 ± 0.029 it 0.007. If the results reported in second reference 
of [31] is confirmed, it would also rule out any significant non-perturbative effects that 
predict strange quark contribution to the spin of proton from pion-nucleon sigma term 
(Tttat which in turn alters the value of g\. Such a contribution can be calculated in the 
framework of chiral quark model. In fact, in this framework it is shown that [32] while 
light sea quark polarization is consistent with zero, the AS = —0.051. If it is added to 
our result,it becomes in line with g^ value obtained from hyperon decay. Nevertheless, 
The crucial point in our model is that there is no room for the perturbative sea quark 
polarization, because the sea of the valon is generated entirely from gluon splitting and for 
the massless quark, the helicity is conserved; and yet the model nicely accommodates all 
of the experimental data with acceptable accuracy. Our finding that Ag^ea is consistent 
with zero is in agreement with the HERMES [18] [19] and SMC collaboration [33] data. 
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This can be seen yet in a different way as follows: 

The moments of sea quarks in a proton can be written as 

Ag(n,g2) = ^[AE(n,Q2) _ Au^(n,Q2) _ Ady{n,Q^)] (25) 

where 

Au„(n, = 2AM^^(ra, Q^) AMu (n) 

Ady{n,Q^) = AM^^{n,Q^)^AMn{n) (26) 

p 

AS(n, g2) = AM^(n, Q^) O [2AMc/ (n) + AMd (n)] 

p p 

where AMu^p (n) are moments of 5Gu,d- For n = 1 they are given by 

p 

r dy6Gu {y) = 0.403, dySGoiy) = -0.409 (27) 

Jo p Jo p 

Therefore, EQ.(25)becomes 

Aq{n,Q^) = (AM^(n,Q2)_AM^^(n,Q2))(2AMc/(n,Q2) + AMD(n,Q2)) (28) 

p p 

and we see that for n = 1 the range of variation for A^(1,Q^) is — 0.016, since 
AM'5(1,Q2) = 1 and AM^^{n = 1) varies between 1 and 0.80 for = [0.283, 10^] 
GeV^; that is, the contribution of the sea quark to the spin of proton is consistent with 
zero. Reference [15] introduces two different polarized valon distribution in an attempt 
to avoid dealing with Aq{l,Q'^) « 0. Such an scheme is quite counter intuitive, because 
it means that singlet and non-singlet quark distribution inside a valon alters the valon 
distribution in a hadron. A notion that is hard to understand within the context of the 
valon model. 

Of course a 2% deviation of our results from the value of g"^ can be attributed to the 
presence of sea polarization which then ought to be generated non-perturbatively . This 
point is not considered here nor have we attempted to fit the data in order to extract 
possible sea quark contribution, in particular that of strange sea. Another possible source 

for this deviation can be a poor determination of 5Guj3_. We also note that a plausible 

p 
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alternative to the full SU{3)f symmetry is a "valence" scenario where SU{3)f symmetry 
is maximally broken which is based on the assumption that the flavor-changing hyperon 
/3-decay data fix only the total helicities of valence quark at some appropriately chosen 
input scale = Qq [17]. 

Our prediction for AE, the total quark contribution to the spin of proton, lies in the 
range of 0.410 - 0.420 for = [2,10]. The variation of AE is due to (marginal) 
dependence of Aq^ in the Next-to-Leading Order; because APj^g ^ 0. This result is also 
compatible with the experimental data of Refs. [17, 18], where for the measure range, 
0.023 <x< 0.6, they obtained AE = 0.347 ±0.024 ±0.066. Further experimental support 
for our findings comes from recently published data from COMPASS collaboration [34] at 
Q2 = 10 GeV^ and 0.006 < x < 0.7. It shows that Auy + Ady = 0.40 ± 0.07 ± 0.05 , 
Au + Ad = 0.0 ± 0.04 ± 0.03, and As ± As = -0.08 ± 0.01 ± 0.03. This result is consis- 
tent with Au = Ad. The conclusion is that: if we accept the validity of both HERMES 
[31] and COMPASS [34] data, then the role of sea quark polarization in calculation of 
polarized structure functions, g'^''^''^ is marginal. We have arrived at the same conclusion 
directly by considering only the QCD processes up to the Next-to-Leading order in the 
valon representation. 

4.1 Role of Gluon, Orbital angular momentum, and the spin of proton 

The phenomenological model that is described here is able to account for the experimental 
data with a good accuracy. However, it still remains to accommodate the spin of proton. 
The spin of valon in the absence of gluon is accounted for by the total spin contribution 
of quarks, as we saw in section 2. The large gluon polarization in valon, however requires 
a sizable negative orbital angular momentum to compensate for the gluon contribution 
beyond Qq. The implication for proton is that Ag rises as increases; being around 0.4 
at = 2 GeV"^ and reaches to 0.7 at = 14 GeV'^. For the same range of Q^, then 
Lz varies from —0.1 to —0.4. These results are plotted in Figure 9. The role of Orbital 
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Figure 9: Gluon, Ag, and orbital angular momentum , L^, components in proton as a function of . 

angular momentum, Lz in a valon is to cancel out the gluon polarization completely, 
but this cancelation in proton is partial. Therefore, it is reasonable to speculate that 
about 60% of the net spin of proton comes from gluon. This is comparable with ~ 50% 
momentum contribution of gluon to the total momentum of proton. We do share the 
opinion expressed in [31] based on the experimental data from HERMES that the quark 
helicities contribute a substantial fraction to the nucleon helicity, but there is still need 
for a major contribution from gluon/orbital angular momentum. We have concluded that 
the orbital angular momentum is needed and arises from the structure of valon in order 
to compensate for the growing gluon helicity and produce a spin-i valon. It also can 
be argued that the valons of nucleon themselves might have a relative orbital angular 
momentum other than zero. Such a situation would amount to the assumption that the 
nucleon is not in the S- state of orbital angular momentum, the case that we have not 
considered in this analysis. In the analysis of section 2 we stated that the magnitude of 
gluon helicity, Ag, depends on the initial scale, Qg, chosen; and hence, so does the values 
of Lz- We believe that the mathematical boundary condition of the model provides a 
reasonable guideline for the choice of Qq- 
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5 Remarks and conclusions 

We have calculated the polarized structure of a valon form QCD processes in the next- 
to-leading order framework. While the valence quark completely accounts for the spin of 
a valon, the presence of large gluon polarization in the valon makes it more complicated, 
requiring a sizable and negative orbital angular momentum. Our finding indicates that 
the sea parton polarization in the valon remains very small, and hence its contribution to 
the spin structure of proton is consistent with zero. This finding is in agreement with the 
experimental results [18], [19]. The picture presented here is capable of reproducing all 
available data on ^^'"'''^(Q^) with good accuracy. We have further calculated the orbital 
angular momentum contribution, and its evolution, to the spin content of proton and 
valon. It appears that the size of gluon contribution to the spin content of proton is 
around 60%, somewhat similar to the momentum contribution of gluon to the momentum 
of proton in the unpolarized case. This value is also sensitive to the initial scale, Qq- 
The model presented here does not have any free parameter, it is free of data fitting and 
solely relies on QCD processes, except the use of phenomenological concept of the valon 
model. Finally, we stress that there is the issue of initial input densities at scale, Qq. 
While the most theoretical analysis and global fits begin with Qq> 1 GeV'^, we have used 
the mathematical boundary conditions of the model and have shown that the results are 
compatible with the experimental data. 
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